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PACS 68 . 08 . -p - Liquid-Solid interfaces 
PACS 62 . 20 . mm - Fracture 

PACS 68 . 35 . Ja - Surface and interface dynamics and vibration 

Abstract. - A tenuous monolayer of hydrophobic particles at the air- water interface often forms 
a scum or raft. When such a monolayer is disturbed by the localized introduction of a surfactant 
droplet, a radially divergent surfactant shock front emanates from the surfactant origin and packs 
the particles into a jammed, compact, annular band with a packing fraction that saturates at a 
peak packing fraction <j>*. As the resulting two-dimensional, disordered elastic band grows with 
time and is driven radially outwards by the surfactant, it fractures to form periodic triangular 
cracks with robust geometrical features. We find the number of cracks N and the compaction band 
radius R* at fracture onset vary monotonically with the initial packing fraction (<f>i n u)- However, 
its width W* is constant for all (fiinit- A simple geometric theory that treats the compaction 
band as an elastic annulus, and accounts for mass conservation allows us to deduce that N ~ 
2irR*/W* — 4:-K(j)RCP /4>init, a result we verify both experimentally and numerically. We show the 
essential ingredients for this phenomenon are an initially low enough particulate packing fraction 
that allows surfactant driven advection to cause passive jamming and eventual fracture of the 
hydrophobic particulate interface. 



X 



The behavior of hydrophobic particles at interfaces 
presents interesting phenomena of fundamental signifi- 
cance at the intersection of interfacial physics, chemistry, 
and continuum mechanics [T][?j , as well as being relevant in 
a variety of applications as particles at interfaces can sta- 
bilize drops and emulsions via jamming [3J. Recent experi- 
ments [US] have demonstrated the two-dimensional elastic 
properties of jammed hydrophobic particulate monolayers 
at interfaces and aspects of their deformation via buck- 
ling and random cracking. These particulate rafts also 
afford a nice system to study dynamics of the formation 
of a jammed solid, and its failure via fracture and flow, 
both phenomena that are hard to analyze in bulk gran- 
ular media. Here wc study both these phenomena at an 
interface using a combination of experiments and com- 
putations and show how a tenuous particulate monolayer 
driven by Marangoni stresses induced by localized surfac- 
tant introduction |16j leads to formation of a jammed solid 
and its eventual failure via a regular cracking pattern. 

Fig. shows a schematic of the experimental setup. 
A clean glass petri dish (diameter 0.14 m) filled with 
distilled water to a height of 0.01 m is placed atop a 



light tablet. Teflon coated hollow glass particles (diameter 
d = 50±10/im, specific gravity 0.25) are introduced at the 
air- water interface to form a particulate monolayer with an 
initial areal packing fraction (f>i n it (defined as ratio of total 
initial particulate area to total interfacial area) that varies 
in the range 0.1±0.01 < (f> init < 0.64±0.01. Owing to the 
protocol followed for particle introduction, (f>i n it cannot be 
controlled, but can be measured (please see Methods sec- 
tion). When a clean steel needle wetted with oleic acid is 
dipped into the water surface at the dish center at time 
t = s, the spreading surfactant pushes the hydrophobic 
particles radially outwards and packs them along an annu- 
lus around the growing particle-free hole as shown in Fig. 
[Tp and [T(d] . The compaction dynamics are imaged with a 
high speed digital camera (Phantom v5) at 600 frames per 
second and analyzed to measure the packing fraction (see 
Methods section) . This allows us to follow the compaction 
band's evolution in terms of the azimuthally averaged ra- 
dial packing fraction 4>e{r, t) = ^ (f>(r, 6, t) d0. 

Immediately following surfactant introduction, the par- 
ticles move slower relative to the surfactant front due to 
subphase drag. After a short time, the particle dynamics 
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Fig. 1: (color online) a) Setup: A clean petri dish with dis- 
tilled water is placed on a light tablet. A camera suspended 
vertically above records compaction dynamics of particle rafts 
at the air-water interface. A steel needle introduces surfactant 
at the interface and leads to a radially propagating surfactant 
shock, b) Snapshot of particle covered surface after surfactant 
introduction shows an annular compaction band around nearly 
regular triangular cracks. The compaction band's leading edge 
is termed the compaction front, while the surfactant shock is 
slightly ahead of the trailing edge. The angle a made by the 
crack face with the radial line from surfactant tip (red lines) 
varies within a narrow range across all cracks for all <j)imt- 



settles into a self-similar form that persists for a while. 
Fig. [2k shows (j>g(r,t), normalized by the peak saturation 
packing fraction <j>* at various times for a representative 
experiment. The compaction band moves radially out- 
wards from the location where the surfactant is introduced 
(r = 0) as time progresses [T7], with the surfactant shock 
front slightly ahead of the trailing edge of the band Rt (see 
Fig. [2] a). As the particles ahead are swept up to form the 
compaction band, the peak packing fraction rises through 
a short transient and saturates at 4>g(r,t)/(f>* ~ 1.0 (see 
Fig. (5^). Thus, the initially tenuous low density interfa- 
cial raft gets packed by the surfactant shock and forms a 
jammed disordered solid when the packing fraction satu- 
rates at 0* (at a time £*) whose value was experimentally 
determined to be slightly below Random Close Packed 
Density (4> RC p = 0.84 in 2D [5]). In Fig. [2Ji, we show 
the positions of the compaction band's leading edge Rl, 
trailing edge Rt, and its width W as a function of time, all 
three of which grow with time as r 3 / 4 (Fig. [5Jd), a scaling 



that persists over almost two decades in time. 

In Fig. [2t>, where one observes Rl, Rt, and W do not 
achieve asymptotic scaling until t ~ 0.4i* . This is also sup- 
ported by cf>g(r,t) I '(f)* evolution [17] which does not start 
rising towards peak saturation packing fraction </>* until 
t ~ 0.4t*. The particles become non-inertial following this 
short transient after which they arc passively advected by 
the surfactant front. However, this subtle initial effect has 
no bearing upon the primary experimental quantities of 
interest which are either static {4>mit), or if dynamical, 
only become relevant at t = t* as discussed below. 

The £ 3 / 4 scaling is consistent with the classical result 
[7] that surfactants from a constant source spread in a 
self-similar form in deep fluid layers i.e. the thickness 
of the viscous boundary layer in the fluid bulk is much 
smaller than the depth of the fluid layer over the dura- 
tion of the experiment. This requires the ratio r/T <C 1, 
where r ~ 0.5 s is the total experimental duration and 
T = H 2 /v = 100 s (H = 10- 2 m is fluid layer depth and 
v = 1 x 10~ 6 m 2 /s is kinematic viscosity of water) is the 
time required for the Blasius boundary layer to span the 
entire depth of the bulk fluid. The ratio r/T - 5 x 10~ 3 
for our experiments, placing them in the deep fluid layer 
regime. The position R s of such a surfactant front for 
uncontaminated surfaces (please see Methods section for 
further details) follows the relation 
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where (water— air) T '(oleic — water) '(oleic— air) 

23.58 xlO -3 N/m is the Harkins spreading coefficient at 
the line of three-phase contact, p = 10~ 3 Pa-s is the dy- 
namic viscosity of the underlying fluid, p = 10 3 Kg/m 3 
its density, and K a numerical coefficient in the range 
0.665 < K < 1.52 [5]. On superposing this predicted scal- 
ing for Eq. Q] with K = 0.94 (solid black line in Fig. ^) 
on the experimental scaling for Rl we see the particu- 
late band propagates like the surfactant shock (Thoreau- 
Rcynolds ridge [5]), and implies the particles behave as 
non-inertial tracers advected by the surfactant flow con- 
sistent with earlier experiments [12] ■ The robust t 3//4 scal- 
ing observed even after <f>g (r, t) / <p* saturates shows the as- 
sumption of non-inertial particle dynamics is a good one. 

As particles are swept up into the compaction band they 
form a jammed disordered solid when (pg(r,t)/<p* — > 1 
saturates at a critical time t*. As the compacted band 
moves radially outward, it fractures to form a periodic 
saw-tooth pattern. There is no observable time difference 
between jamming and fracture onset, which is consistent 
with the fact that the elastic strain before fracture in the 
jammed solid of nearly rigid grains is likely to be very 
small. Fracture thus starts at t = t*, and we therefore de- 
fine R* = Rrit*) and W* = W(t*) as the inner radius and 
width of the compaction annulus respectively at fracture 
onset. 

Beyond the time t > t* , we almost never see the for- 
mation of any new cracks so that the number of cracks 
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Fig. 2: (color online) Experiment: a) Azimuthally averaged ra- 
dial packing fraction <j>e(r,t)/<f>* vs. r shows compaction shock 
evolution at t = 0, 2, 2.3, 2.67, 3, 3.3, 4.17, 8.3, 16.7, and 25 
x 1CP 2 s from left to right. Vertical dashed lines (red) mark 
the compaction band's shock front (Rl), rear {Rt), and width 
(W) at t = 0.083 s. t* (indicated at top left) represents the in- 
stant when peak packing fraction saturates at </>*. b) Position 
of Rl (black circle), trailing edge Rt (red square) and width 
W (blue plus) vs. time t in log-log scale. The solid (black) 
line for R s scaling (Eq. [T]for K — 0.94) exactly superimposes 
Rl scaling, c) Measured no. of cracks N (solid black circles), 
Eq. [5] (solid red squares), and Eq. [3]vs. <j>i„it- Inset: W* /d is 
constant for all cbinit- 



N formed remains constant in a given experiment, even 
though the cracks grow dynamically. We find with increas- 
ing initial particulate packing fraction (j>i n it, the number 
of cracks decreases monotonically (see Fig. |2fc) as does the 
critical radius R* . This is because the compaction annulus 
forms and jams at an earlier time and smaller radius with 
increasing <pi n u ■ However the compaction band width W* 
exhibits no dependence on <\>initi but as we explain later, 
it does depend on the particle diameter (see Fig. [3]). Since 
cracks in the compaction band relieve strains over a scale 
comparable to the band width W* , we expect that the 
number of cracks 



N ~ 2irR*/W* 



(2) 



Additionally, assuming the initial particle distribution is 
uniform and the jammed solid is random close packed, 
mass conversation dictates that the particulate area within 
the compaction annulus at t* equals the particulate 
area within a circular radius (R* + W*) at t = 0, so 
((Prcp7t[(R* + W*f-{R*) 2 ] = <t> ini MR* + W*) 2 ) so that 
(f>RCp/(f>init = R*/2W*, thus yielding 



N ~ 2nR*/W* ~ ^RCp/^inU 



(3) 



In Fig. [2J;, we plot the experimentally measured values for 
N vs. ^C 2 4>BCp/<t>init, with C 2 = 0.7. 

The agreement of Eq. [5] with experiments suggests 
that the continuum description holds well for this gran- 
ular system, a fact also supported by Fig. [5£ inset where 
W*/d ~ 50 ± 10 at all 4>i n it- However, Eq. [3]dcrivcd from 
mass conservation arguments shows only partial agree- 
ment. We emphasize that this relation is based on sev- 
eral idealized assumptions, and holds only for intermediate 
values of <\>imt- Firstly, we drop the quadratic term dur- 
ing expansion of the mass conservation relation leading to 
Eq. [3] under the assumption that W* <C R* (or alter- 
natively N 3> 1), which holds only when 4>i n u < 4>RCP- 
Secondly, we assume an initially uniform particle distri- 
bution. In reality, the meniscus formed by water with the 
petri dish wall repels hydrophobic particles towards the 
center. This is clearly observed as a drift in <j)g(r,t) for 
r > Rl in Fig. [2^,. Finally, we assume the compaction 
annulus jams at <Pb,cp- In reality, frictional or attractive 
inter-particle interactions can stabilize a granular pack be- 
low 4>rcp- Whereas the teflon coat on particles ensures 
minimal inter-particulate friction, we do observe particle 
clustering suggesting attractive interactions are at play. 

We expect both Eq. [5] and [3] will fail in the limiting 
cases as (f>i n u — > 0,4>RCP- in th° dilute limit, Eq. [5] and 
E] suggest N — > oo, which cannot be true. Instead, A" is 
strongly influenced by disorder in the packing, since the 
crack size is only a few grain diameters, and furthermore 
A" < 2wR* jd because a crack cannot be smaller than par- 
ticle diameter d. The continuum description however fails 
well before this limit is reached. In the opposite limit 
as 4>init — > 4>RCP, the particulate layer is either already 
jammed, or does not have to be packed much before it 
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• Teflon coated hollow glass spheres 
■ Pliolite VTAC-L particles 
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Fig. 3: (color online) Dimensionless compaction width W* /d 
vs. particle diameter d for d = 50±10/mi (Teflon coated hollow 
glass spheres), 90±20Mm, 130±2(kt m, 180±30/xm, 230±2(ktm, 
330 ± 30/mi, and 460 ± 30/xm (all for Pliolite) shows W* scales 
with particle diameter. Horizontal error bars represent parti- 
cle dispersity, whereas vertical error bars represent variability 
in measured W* over 20 experimental runs for Teflon coated 
hollow-glass spheres, and 4 experimental runs each for Pliolite 
particles. 



jams into a solid, so that t*,R*, and W* are all poorly 
denned. Also, since particles are already constrained and 
no free space is available for cracks to open up and expose 
surfactant to the air, cracks must proceed through local 
re-arrangement of particles which leads to branching and 
kinking instabilities [4]. 

As shown in the inset of fig. [5J;, we found the di- 
mensionless compaction width W*/d is independent of 
the initial packing fraction (pmu- This critical width W* 
depends upon the particle diameter d, and the surface 
tension contrast ratio A7/7 (dimensionless Marangoni 
stress) . Here 7 is the surface tension of the oleic-air inter- 
face (7 = 32.8 x 1CP 3 N/m). Since the critical strain 
for fracture must depend on these parameters, dimen- 
sional analysis and St. Vcnant's principle suggest that 
e = (W* /d)(A r y/j). Additional experiments with varying 
particle diameter were performed to verify this. In addi- 
tion to teflon coated hollow glass spheres, we also used 
polydisperse Pliolite hydrophobic particles of various 
diameters. Fig. [3] shows the dimensionless compaction 
width W* jd vs. particle diameter d has zero slope indi- 
cating W* scales linearly with particle diameter d. 

Having considered the formation and number of cracks, 
we now analyze the crack geometry. Since all dimensions 
of the compaction annulus exhibit self-similar scaling (see 
Fig. [5b), the particle- free area opened by a crack too must 
scale self-similarly. Thus we expect the crack growth to 
also be self-similar, since they are driven by stresses in 
this growing annulus; the toothed crack front with straight 
edges naturally fits these constraints. An individual crack 
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Fig. 4: Molecular dynamics simulation for formation and fail- 
ure of a jammed solid. A compaction band is shown for (a) 
A7 = 3k /d and (b) A7 = 6k /d. The non-dimensional sys- 
tem radius is r/d = 400 (in experiment the petri dish radius 
r/d — 1400). (c) Simulated profiles of azimuthally averaged 
packing fraction </>e (r, i) show compaction shock evolution with 
increasing Rl- (d) Simulated number of cracks is compared to 
Eq. [3] vs. 4>init- In (c) and (d) A7 = 3k /d. 



is thus all that we need characterize, and this can be done 
in terms of the angle (a) between the crack face and the 
radial direction from the point of surfactant introduction 
(Fig. Q~p) for all cracks and all values of (f>i n it- The dis- 
tribution of a is quite sharp with a coefficient of variation 
of 34.0°/6.97° = 4.9, that is relatively constant in time as 
cracks grow with the diverging compaction band. 

We now turn to a quantitative description of the dynam- 
ics of compaction, jamming and fracture to complement 
the qualitative description of the phenomenon at a seal- 
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ing level. We model the initially tenuous raft as a planar 
system of hard particles with a pair-wise attractive force 
F = k at separation r < d/10, and F = for r > d/10 
with the particle diameter assumed to be uniformly dis- 
tributed in the interval [0.8c?, 1.2c?]; the actual interaction 
potential [TU] is more complicated, but the simple rep- 
resentation described above captures all the qualitative 
trends. @ The initial configuration is created by first plac- 
ing the particles randomly and then relaxing the system 
which leads to some reordering initially. A symplectic Eu- 
ler scheme built into a molecular dynamics simulation |13j 
is then used to solve for the damped Newtonian dynamics 
of the particles. The role of fluid drag on the particles 
is complicated by the presence of a boundary layer. To 
mimic this accurately with an implicit fluid, we assume an 
outward radial flow with velocity U r = inside a radius 
R s given by eq. (|T|) and zero elsewhere [5], that simulates 
the spreading of the surfactant and the resulting fluid flow. 
We set R s {i) = d initially, where i = [d/AT^) 1 / 4 ] 4 / 3 fol- 
lowing cq. ((TJ), to avoid a divergence of velocity at t = 0; 
our simulations therefore start at t = t (the long time dy- 
namics are independent of t). Motion of a particle with 
a velocity v relative to that of fluid is opposed by a drag 
fn = ftdv, where we introduce an effective viscosity: 

v = \—i — r d - ( 4 ) 

y <Pinitt 

which follows from the requirement that a radial inte- 
grated pressure difference across the compaction band 
must be of the order A7. It can be derived by using a 
similar boundary layer argument that has been used to 
derive eq. (QJ itself [8]. According to the Blasius bound- 
ary layer theory, drag exerted per unit width of a flat 
plate is D w U 3 ^ 2 y/lfip where U is the flow velocity and I 
width of the plate in the direction of the flow [TT]. Here 

1 = w « Wsf and u ~ y icldin g D ~ A ^vS- 

On the other hand, number of particles per unit width 
of the band is approximated by Multiplying this by 
fit and assuming all the particles moving with velocity U 

leads to total drag D v ~ ^,??^'"' t exerted on the raft. By 

o P dttpncp ■' 

requiring D p = D one obtains eq. 

In the non-inertial regime explored here, the dynamics 
of the model depends on the radial pressure relative to 
attraction between the particles, which can be controlled 
by A7 and k. For A7 = 3k /d the model displays forma- 
tion of the compaction band and its fracture in a manner 
similar to the experiments (Fig. 0^, [IB])- In the model 
we can accurately measure that the compacted packing 
fraction saturates close to c/>rcp- Fig- Ht US] shows time 
evolution of the azimuthally averaged packing fraction in 
good agreement with experiments (Fig. The scaling 
of N vs. 4>i n it (Fig. 0J1) from simulations also compares 

2 We also tested harmonic interaction potential without notable 
difference in the results. 



well with experiments (Fig. as do the deviations ex- 
pected when 4>init —> and 4>i n it — > 4>RCP a s discussed 
earlier. In addition, the N depends on A7 through the 
pre- factor C2 (see Fig. |2t and|4ji). The model suggests 
that the crack depth is controlled by the radial pressure 
such that for increasing A7 the crack depth decreases, see 
Figs. @Ji,b. Our simulations also allow us to measure the 
crack angle (Fig. 0} a = 30.5° ± 3.5° for A7 = 3k/d and 
a = 33.9° ± 7.2° for A7 = 6k/d, results which compare 
well with experiments. 

The simple molecular dynamics scheme adopted here 
shows the basic features of the experimental phenomenon 
can be captured with a simple attractive interaction be- 
tween particles, irrespective of its actual form. However, 
knowing the actual form of the interaction will better al- 
low one to understand the individual and collective par- 
ticle dynamics. Additionally, understanding the role of 
rough contact lines and particle anisotropy also requires 
careful future study. 

Our experiments and simulations have allowed us to un- 
derstand the formation and failure of the resulting com- 
pact, disordered solid in terms of a structural parameter, 
the initial packing fraction 4>i n it, driven by a differential 
surface tension A7. Since the two features are common to 
many problems involving the mechanics of disordered ma- 
terials, this system might serve as a paradigm for further 
studies in amorphous solids. 

Methods. — Preparation: Standard cleaning proce- 
dures |14j were followed to ensure an impurity free setup. 
The pctri dish was washed in dilute Sulphuric acid, rinsed 
in distilled water, then baked dry at 100°C for 30 min- 
utes, followed by exposure to Ultraviolet radiation in an 
oxygenated environment to break up residual organic im- 
purities. The needle used for surfactant introduction was 
washed in Ethanol, rinsed in distilled water and flame 
treated prior to each experiment. All experiments re- 
ported here were performed with pure Oleic acid (no dilu- 
tion with an organic solvent). The amount of surfactant 
introduced had no bearing upon the results. 

Hollow glass microspheres composed of borosilicate 
glass d were coated with a thin layer of Polytetrafluo- 
roethylene (Teflon) via molecular vapor deposition. Prior 
to the experiment, the particles were washed in Ethanol 
and rinsed with distilled water and baked dry at 100°C. 
Particles were introduced by puffing them in air and allow- 
ing them to naturally settle onto the interface. The initial 
packing fraction (j>i n u cannot be controlled in this particle 
deposition scheme. A number of experimental runs were 
performed by varying the approximate number of particles 
released which indirectly permitted us to span a range of 

Surfactant: The validity of Eq. [TJrcquires that the prop- 
agating surfactant front be introduced from a source of 
constant concentration. This requirement was experimen- 

3 Trelleborg Emerson & Cumming Eccospheres. Product No. W- 
25, mean size 50/^m, density 0.25 g/cc, flotation 95% bulk vol. 
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tally met by introducing the surfactant from a point source 
(needle) with a concentration well in excess of the surfac- 
tant's Critical Micelle Concentration (CMC), which also 
ensured the surfactant surface tension 7 remained con- 
stant across all experimental runs and was not a function 
of surfactant concentration as also confirmed in 0]. For 
an Oleic acid molecule with cross-sectional area of 20 A 2 , 
the total number of molecules required to form a mono- 
molecular layer at the petri dish's air-water interface is 
7.7 x 10 molecules which translates to an approximate 
total Oleic acid mass of 0.36 fig (at Oleic acid density of 
282.46 g/mole). The average droplet mass introduced by 
the needle was 7.66 mg (~ 8.5 /zl volume per droplet), 
more than 2 x 10 4 times the quantity required to form a 
monolayer. Hence we are confident the surfactant spreads 
as a thick surfactant layer and its surface tension is con- 
centration independent. We also conducted hve indepen- 
dent measurements each for the surface tension at the air- 
water, air-oleic acid, and water-oleic acid film (formed by 
introducing droplet with the needle) interfaces using the 
Wilhelmy plate method which confirmed the surface ten- 
sion for surfactant film was concentration independent. 

Image Analysis: The high speed camera (Phantom v5.0 
camera, exposure time: 150 (is with 28 mm Nikkor man- 
ual focus wide angle lens, aperture setting: f/5.6) recorded 
light transmitted through the particulate layer providing 
high contrast images of dark particles in a bright back- 
ground. All images were collected under same illumi- 
nation conditions i.e. DC illumination source intensity, 
camera lens aperture, and exposure time were kept con- 
stant across all runs. A digital snapshot of the background 
(petri dish with distilled water prior to particle introduc- 
tion) was subtracted from images with particle dynam- 
ics thereby removing inevitable minor spatial illumination 
inhomogeneities. Image analysis algorithms for particu- 
late area measurement were developed in-house and first 
calibrated against particles of known area (measured un- 
der a microscope) that were introduced at the interface. 
Tests performed against images obtained from molecular 
dynamics simulations provided the error bars we present 
in Fig. [2J The measurement error for (j)init was estimated 
at S<pi n it = 0.01. For dynamical measurements, control 
tests against molecular dynamics images yielded a lin- 
ear increase in error for <j> > 0.7 reaching a peak error 
8(j>e{r,t) = 0.04. Given this higher error at high pack- 
ing fractions, we are unable to confirm whether the com- 
paction band saturates at 4>rcp = 0.84 or at a lower value 
due to attractive interactions as discussed earlier. In any 
event, this has no bearing upon the results since we only 
seek to learn the instant t* when the peak packing fraction 
saturates to a maximum value <jf heralding the formation 
of a jammed solid. 
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